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We generalize simplicial minisuperspace models associated with restricting the topology 
^ ■ of the universe to be that of a cone over a closed connected combinatorial 3— manifold by 
considering the presence of a massive scalar field. By restricting all the interior edge lengths 
and all the boundary edge lengths to be equivalent and the scalar field to be homogenous on 
the 3— space, we obtain a family of two dimensional models that include some of the most 
relevant triangulations of the spatial universe. After studying the analytic properties of the 
action in the space of complex edge lengths we determine its classical extrema. We then 
obtain steepest descents contours of constant imaginary action passing through Lorentzian 
classical geometries yielding a convergent wavefunction of the universe, dominated by the 
contributions coming from these extrema. By considering these contours we justify semiclas- 
sical approximations based on those classical solutions, clearly predicting classical spacetime 
in the late universe. These wavefunctions are then evaluated numerically. For all of the 
models examined we find wavefunctions predicting Lorentzian oscillatory behaviour in the 
late universe. 
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1 Introduction 



The quantization of gravity, (QG), is perhaps the most important problem in theoretical 
physics today. Among the several avenues proposed to achieve that goal, one of the most 
productive has been the sum over histories formulation. In such a formulation an amplitude 
for a certain state of the universe is constructed by summing over a certain class of physically 
distinct histories that satisfy appropriate boundary conditions, weighted by their respective 
action. Problems with the convergence of the path integral for gravity make it convenient 
to use the Euclidean version of the sum over histories. 

In the (Euclideanized) classical theory of gravity, i.e., general relativity, (GR), a classical 
history is a Riemannian manifold, that is, a smooth manifold (M", A) with smooth structure, 
(atlas). A, endowed with a Riemannian metric g. Topological manifolds M", with their 
local homoeomorphisms to R"^, are the mathematical implementation of GR's Equivalence 
Principle. The smooth and metric structures are essential for the definition of the most basic 
physical concepts like distance, curvature, field derivatives, etc. 

When building a quantum version of GR there are several reasons why we should consider 
the generalisation of this concept of history. Any sum over histories can only be implemented 
once we specify what are the elements of the set of (physically distinct) histories in QG we 
have chosen to consider. Such choice is constrained by several criteria that any reasonable 
history in QG must meet. The most fundamental are: 

• It should be based on a finite dimensional topological space X" endowed with a smooth 
structure A. 

• The above mentioned space must be metrizable. 

• The set of physically distinct histories {X''\ A, g^^^} considered in the sum must be 
algorithmically decidable, AD, and classifiable, AC. 

• The action functional used should coincide with the usual Einstein action for manifold- 
based histories. 

Given such a set of histories it is then possible to write a probability amplitude for 
topology change between a set of boundary (n — 1)— manifolds dX"^~^ in the same cobordism 
class, as 



G[dX--\h]= / /^^7^.e-^[^"'^'^-l (1.1) 

(X",A) 

where the sum is over all physically distinct histories, {X^^A^g^y}, that have the ap- 
propriate boundaries ((9X"~^), and induce the desired {n — 1) metric, on those same 
boundaries, and / is the Euclidean action associated with each history. 

Any successful theory of QG must have GR as its low energy limit, and so its space of 
histories must contain the space of classical, manifold-based histories. However, although 
necessary there are several indications that they are not sufficient. To start with, the set 
of n— dimensional smooth manifolds is not AC for n > 4, and there is no known algorithm 
to decide whether or not a generic topological space is a manifold in 4 dimensions. It was 
proven there is no way to do so for higher dimensions. 

Furthermore, there are cases in which formal descriptions of histories based on the clas- 
sical configurations of a theory, do not necessarily correspond to the precise mathematical 
definition of the space of histories needed in order to make Euclidean sums over histories 
both well defined and yielding the correct quantum mechanics, |jl[. On the other hand these 
histories cannot be based on very pathological topological spaces, because it would be im- 
possible to concretely define the action associated with those generalized topological spaces. 
This requires the definition of concepts like distance, volume and scalar curvature. In order 
to define distance, the topological spaces must be metrizable, and it can be shown that a 
notion of curvature can be introduced in any metrizable space where such notion of distance 
exists. 

However, such spaces must also have uniform dimension, otherwise it would not be clear 
how to weight their contributions in a sum over histories, since the form and properties of 
the action depend on the dimension of the space. As pointed out by Scheleich and Witt in 
0, these restrictions are still not enough to eliminate all unsatisfactory topological spaces, 
and they go on to discuss in some detail what other restrictions are necessary. 

A good candidate for this generalized space of histories as been suggested in 0. They 
propose generalizing from manifold-based histories to conifold-based histories. Briefiy a coni- 
fold is a topological space that is like a manifold everywhere except in a discrete, countable 
set of points S, (called singular points), that do not have any neighbourhood homeomorphic 
to an open ball of BJ^, but to a cone over some closed connected (n — 1)— manifold (other 
than S'""^, of course). 



In one and two dimensions, manifolds and conifolds coincide. But for higher dimen- 
sions the set of conifolds more than contains the set of manifolds. The fact that the set 
of non- manifold like points, S, in any conifold is discrete and countable allows the direct 
generalization of all basic geometrical concepts, (such as geodesies, curvature, etc.) from 
manifolds to conifolds by continuation. 

The problem now is how to concretely implement such a sum. To do it we must not only 
define the space of histories we will be considering but also obtain a finite representation 
for them. A simplicial formulation of QG in terms of Regge calculus aims to provide one 
such representation. Following |^ by defining combinatorial manifolds and conifolds in the 
appropriate way, we can be sure that in less than seven dimensions there is a one to one 
relationship between smooth manifolds/conifolds and their combinatorial counterparts. 

These combinatorial counterparts are thus a finite representation of the smooth topo- 
logical spaces they are associated with. Furthermore this one to one relationship allows us 



to substitute the sum ( |1 . 1| ) for a sum over those simplicial complexes, which is easier to 
concretely implement. 

However, the computation of the full sum is not at present viable, and so we restrict 
ourselves to a simplicial minisuperspace approximation for the wavefunction of the universe 
for some simple but fairly general models. Such calculations were initiated in ||^, Q and 
extended in ^ and We aim to further extend these models and verify if they yield a 
wavefunction predicting classical spacetime. 

Our simplicial minisuperspace approximation consists in restricting the topology of the 
4 — D universe to be that of a simplicial cone over a 3 — D closed connected combinatorial 
manifold, which is to be seen as a triangulation of the spatial 3 — D universe. Furthermore, 
its geometry must be such that there is a single type of boundary edge ,with square length Sb, 
and internal edge with squared length, Sj. We shall also consider that there is a massive scalar 
field present, but that it takes the same value in all the vertices of the spatial 3 — D universe, 
which is analogous to the requirement of the scalar field to be homogeneous (p = 4'{t) in the 
usual continuum minisuperspace models, |TD . 



As pointed out in [Q, we shall see that the only way of obtaining a wavefunction pre- 
dicting a late universe exhibiting classical spacetime like our own, (i.e. a wavefunction with 
Lorentzian oscillatory behaviour), is to consider contours of integration over complex valued 
interior edge lengths, Sj. This forces us to study the minisuperspace action as a function of 
complex variable, which is rather involved because of the multivaluedness of the action. 



We then try to find convergent integration contours yielding appropriate wavefunctions, 
i.e., that predict Lorcntzian classical spacetime in the late universe, by allowing a semiclas- 
sical approximation based on Lorentzian classical solutions to be a good approximation of 
the full wavefunction when the universe is large. 

2 Simplicial Quantum Gravity 

2.1 Combinatorial Manifolds and Conifolds 

Before we introduce our simplicial minisuperspace model it is convenient that we review 
some basic definitions of simplicial geometry. 

Definition 2.1 A simplicial complex {K, \ K \) is a topological space \ K \ and a collection 
of simplices K ,such that 

• \ K \ is a closed subset of some finite dimensional Euclidean space. 

• If a is a face of a simplex in K, then a is also contained in K . 

• If aa and are simplices in K , then a a is a face of both a a and (76 . 

• The topological space \ K \ is the union of all simplices in K. 

Definition 2.2 A combinatorial n— manifold Ai^ , is an n— dimensional simplicial complex 
such that 

• It is pure. 

• It is non branching. 

• Any two n— simplices can be connected by a sequence of n— simplices, each intersecting 
along some {n — 1) — simplex. 

• The link of every vertex is a combinatorial {n — 1)— sphere. 

Note that there are simplicial complexes that are homeomorphic to topological manifolds 
but are not combinatorial manifolds. The definition of combinatorial manifold carries more 
structure than simply the topology. A similar definition can be made for combinatorial 
conifolds C", by simply replacing the previous last condition by 



• The link of every vertex of is a closed connected combinatorial {n — 1)— manifold. 



Thus, like in the continuum framework, a general combinatorial conifold is a combinato- 
rial manifold everywhere except possibly at a countable set of vertices. 

Combinatorial manifolds and conifolds are in a similar situation to smooth manifolds 
and conifolds, in the sense that they both have additional structure relative to their more 
general topological counterparts. Such structure allows integration and differentiation to be 
well defined as essential for any physical applications. In order to see the connection between 
the smooth and combinatorial structures we need to introduce one more definition: 

Definition 2.3 A combinatorial triangulation of a manifold, M", consists of a combinato- 
rial manifold TVI", and an homeomorphism t :| TVI" |— M". 

An analogous definition holds for conifolds. 

Following [0] it can be shown that every smooth manifold and conifold admit combi- 
natorial triangulations. On the other hand, it can be also shown that every combinatorial 
manifold admits a piecewise-linear, (PL), structure, i.e. a PL-atlas, {(K, 0a)}aeA, such that 
the mappings 

0,0-1 : (t>a{Va n H) 06(K n H) (2.1) 

are PL mappings between subsets of i?" . Finally it is known that in less than seven dimen- 
sions every manifold with a PL-structure has an unique smoothing. A similar reasoning can 
be applied to conifolds, and so we are able to conclude that 

In less than seven dimensions, smooth manifolds/conifolds, (M"/C"), uniquely corre- 
spond to combinatorial manifolds/conifolds, (Al'^/C"). 

Obviously the topological space underlying each smooth manifold/conifold may have 
several inequivalent triangulations, what this result says is that when we specify a smooth 
structure on that space it will correspond to a unique combinatorial triangulation, i.e., a 
triangulation based on a combinatorial manifold/conifold, and not just any simplicial com- 
plex. This unique correspondence implies that in less than seven dimensions the information 
about the smooth structure of the space is carried by its combinatorial triangulation. So 
we see that the topological part of the sum over histories in ( |1.1|) , can be recast in terms of 



simplicial representatives of the "continuum" spaces. However in order to fully translate an 
heuristic expression such as ( |1 . 1| ) , into a concrete sum over simplicial histories we still need 
to associate a metric and an action to each simplicial complex to be considered. Note that 
up to now we have not specified any kind of metric information associated with simplicial 
complexes. Once we have fixed the topology of the underlying simplicial complex, the most 
convenient way to attach metric information to it, is to use Regge calculus. 

2.2 General Regge Formalism 

A convenient way of defining an n— simplex is to specify the coordinates of its (n + 1) vertices, 
a = [0, 1, 2, ...,n]. By specifying the squared values of the lengths of the edges Sij, we 
fix the simplicial metric on the simplex. 

/ N Soi + Soj — Sij 

9ijisk) = ^ (2.2) 

where i,j = 1, 2, ..n. 

So if we triangulate a smooth manifold M endowed with a metric g^i, by a homeomorphic 
simplicial manifold Ai, the metric information is transferred to the simplicial metric of that 
simplicial complex 

= ^ ~ (2.3) 

In the continuum framework the sum over metrics is implemented through a functional 
integral over the metric components {g^u{x)}. In the simplicial framework the metric degrees 
of freedom are the squared edge lengths, and so the functional integral is replaced by a simple 
multiple integral over the values of the edge lengths. But not all edge lengths have equal 
standing. Only the ones associated with the interior of the simplicial complex get to be 
integrated over. The boundary edge lengths remain after the sum over metrics and become 
the arguments of the wavefunction of the universe. 

J Dg,,{x) j D{s,} = 11 J dfi{si) (2.4) 

In the simplicial framework the fact that the geometry of the complexes is completely 
fixed by the specification of the squared values of all edge lengths, means that all geometrical 
quantities, such as volumes and curvatures, can be expressed completely in terms of those 
edge lengths. Consequently the Regge action (the simplicial analogue of the Einstein action 



for GR) associated with a complex of known topology can be expressed exclusively in terms 
of those edge lengths. 

The Euclideanized Einstein action for a smooth 4— manifold M with boundary 9M, and 
endowed with a 4— metric, g^j_p, and a scalar field $ with mass m, is 



JM IGttG JdM SnG 

I JM 

where K is the extrinsic curvature. 

Its simplicial analogue will be the Regge action for a combinatorial 4— manifold, A^, with 
squared edge lengths {sfc}, and with a scalar field taking values {0^} for each vertex v of 



—2 2A 



cr^ 0-1 = [U] 

2^ - ^v.. -^ '/'j- 



where: 

• (Ta: denotes a fc— simplex belonging to the set of all fc— simplices in M.. 

• 9{al), is the deficit angle associated with the interior 2— simplex = [ijk] 



B{cy\)=2-K- Qd{cy\.^^) (2.5) 

and 0(i{(72, (X4) is the dihedral angle between the 3— simplices (73 = [ijkl] and a'^ = 
[ijkm], of (T4 = [ijklm] that intersect at a^. Its full expression is given by 0. 

• '^/'(crg) is the deficit angle associated with the boundary 2— simplex cTg: 



H4) = ^- E 

cr465t(o-|) 



(2.6) 



• Vk{crk) for /c = 2,3,4 is the A;— volume associated with the fc— simplex, cr^, and once 



• V4((Ji), is the 4— volume in the simplicial complex M., associated with the edge ai, i.e., 
the volume of the space occupied by all points of Ai that are closer to o"i than to any 
other edge of Ai. The same holds for V^lJ) where j represents all vertices of A4. 

It is easy to see that both V4(c"i) and V^i^j), can be expressed exclusively in terms of the 
edge lengths {sk}- All these expressions remain valid if we consider smooth conifolds and 
their combinatorial counterparts. 

So we see that any reasonable history in QG of the type {X^, A, g^,^,^), where rep- 
resents either a topological manifold or conifold endowed with a smooth structure A, metric 
g^i^, and in the presence of matter fields represented by $, has an unique simplicial analogue, 
(A"^, {sk}, {^j})- This allows us to recast the heuristic expression ( |l.l| )in terms of this finite 
representation as: 



• {si} are the squared lengths of the interior edges 

• {sb} are the squared lengths of the boundary edges 

• {(pi} are the values of the field at the interior vertices 

• {(ph} are the values of the field at the boundary vertices 

Although the functional integral over metrics has been written explicitly in terms of the 
edge lengths, this expression is still heuristic because we still need to specify the list of 
suitable simplicial complexes Af^ we intend to sum over, the measure, and the integration 
contour to be used. One way to avoid the problems in defining one such list, namely the 
problems of algorithmic decidability and classifiability, (that have been discussed extensively 
in PI ) is to evaluate the sum approximately by singling out a subfamily of simplicial histories 
described by only a few parameters and carrying out the sum over these histories alone. 

An example of this is to adopt a simplicial minisuperspace approximation. We now 
describe in some detail the minisuperspace model we shall consider. 



again their explicit expressions in terms of the squared edge lengths are given by ||^. 




(2.7) 



where 



3 Simplicial Minisuperspace 



We shall reduce our attention to a significant subfamily of simplicial histories characterized 
by the following restrictions: 

3.1 Topological Restrictions 

We shall consider that the universe is well approximated as a simplicial cone = a * Ai^ 
over a 3— dimensional closed combinatorial manifold Ai^, that is taken to triangulate the 
spatial 3— dimensional universe. Models of this kind have already been considered in [^]. 

The most distinctive feature of such a model is the simplifications introduced by the 
existence of only one interior vertex, the apex of the cone. 

The consequences of this restriction are 

• Note that the simplicial 4— complex, C^, in general will not be a combinatorial 4— mani- 
fold but a combinatorial 4— conifold, because all the vertices in are manifold-like 
points, except the apex, a, whose link is L(a) = and in general Ai^ will not be a 
combinatorial 3— sphere. 

• The fact that Ai^ is closed means that the only boundary of will be the itself: 



So the combinatorial 3— manifold is to be seen as a triangulation of the spatial 
3 — D universe. 

• The cone structure of the reflected in the fact that there is only one interior vertex 
(the apex a) means that all vertices in will be boundary vertices of C^. So if Np{}C^) 
is the number of p— simplices in the complex /C" we see that since = a * Ai^, then 



And in general there will be two kinds of p— simplices (p = 1,2,3), in C^, the ones 
that exist originally in A^^, which will be the boundary p— simplices of C^, and the 
p— simplices generated as cones over the {p — 1)— simplices of Ai"^ with apex a. These 
will be the interior p— simplices of C^. So in general 



(3.1) 



No{C^) = No{M^) + 1 



(3.2) 



and 

7Vo(C4) = NoiM') + 1 
Ni{C^) ^ Ni{M^) + No{M^) 
N^iC*) = N^iM^) + N,{M^) 
N^{C^) ^ N^{M^) + N2{M^) 

• The Euler characteristic of a simphcial complex /C" is 

X(/C") = iVo(/C") - iVi(/C") + iV2(/C") - iV3(K:")... ± N^^ilC^) 

with + if n is even and — if n is odd. On the other hand, since A^^ is a closed 
combinatorial manifold, then its Euler characteristic must vanish, and so we have 

No{M^) - Ni{M^) + N2{M^) - NsiM^) = (3.3) 

• Since is a closed, pure, non branching complex then all its 2— simplices must belong 
to exactly two 3— simphces of 

N2{M^) = 2N3{M^) (3.4) 

In the tabic below we list the values of A^O) and for some of the most relevant closed 
connected combinatorial 3— manifolds 





No 


Ns 


Cm=0 
^crit 




5 


5 


29.31 




10 


30 


5.61 


£(2,1) 


11 


30 


3.54 


£(5,1) 


15 


89 


-0.84 




15 


90 


-0.31 



where 

• a4 is the simplest triangulation of S^, [|]. 

• S"^ X and are simple triangulations of the spaces S"^ * and constructed in 

• £(p, 1) are simple triangulations of their respective Lens spaces, L{p,l) used in |]5|. 
The meaning and relevance of the values in the last column will be explained later. 
3.2 Metric Restrictions 

Up to now we have concentrated on the restrictions on the topology of the simplicial space- 
times that characterize our minisuperspace models. The restrictions on the metric degrees 
of freedom are as important. 

By assuming a cone-like structure = a * A4^, we see that all the interior edges of 
are of the same type ,i.e., one of the vertices is a boundary vertex belonging to and the 
other is the (single) interior vertex of C^, its apex. 

If we label the interior vertex as and the other (boundary) vertices of as 1,2, 
NolAi^). Then the cone-like structure of leads to all interior edges being of the same 
form [0,a], with a = 1,2, .., No{M^). 

So it makes sense to introduce the restriction that all interior edges have equal lengths 
whose squared value is denoted Si = soa- A similar assumption is made with respect to the 
boundary edge lengths, i.e., we consider them all to be equal to a common value Sij = Sb, 
with t,j = l,2,..,iVo(-M'). 

We are well aware that these simplifications greatly reduce the scope of the model, namely 
because they result in there being only one kind of 4— simplex in the complex. Thus any 
real classical solution will necessarily be either purely Lorentzian or Euclidean. However, 



it is well known that in the analogous continuum minisuperspace model, complex classical 
solutions not only exist but in the late Universe even dominate the path integral, |]TU|. In 



order to begin to circumvent this problem we are currently studying a similar model with 
two different interior edge lengths. 

3.3 Scalar Field 

The simplifications assumed in respect to the edge lengths makes it natural to assume that 
the scalar field is spatially homogeneous. So we assume that the scalar field takes the same 
value 0b for all boundary vertices of C^. The value at the interior vertex, 0j, is independent. 

3.4 Minisuperspace Wavefunction 

We can now concretely implement a simplicial minisuperspace approximation to the wave- 
function of the universe of the type (|2.7| ), as 



^[M^Sb,(l)b]= dsMie-^^"*^''''''"'^''^"^ (3.5) 



The Regge action for this minisuperspace can now be calculated. For simplicity we 
introduce rescaled metric variables: 

e = - (3.6) 

Sb 



S = ^ (3.7) 



P 

where = PA/3, and P = IGirG is the Planck length. We shall work in units where 
c = h = l. 

Then the volume of the 4— simplices in is 



The volume of the A^i(A^ ) internal 2— simplices, in C is 



P 



= i^,S^i - 1/4 (3.9) 



The volume of the 2N^{}A ) boundary 2— simplices, in C is 



V2{al) = (3.10) 



The volumes of the internal and boundary 3— simplices of C are, respectively 



3 



V^3(aD = Y^^-^'^Vse - 1 (3.11) 

Vs{4) = ^^^/^ (3.12) 
and so the dihedral angle associated with each internal 2— simplex is 

eiat.) = arccos^|— ^; (3.13) 
for the boundary 2— simplices we have 

9{a',) ^ arccos ^^/g^ - (3.14) 

With respect to the matter terms, the kinetic term vanishes when the edges o"2 are 
boundary edges. The only non vanishing contribution comes from the internal edges a2 — 

m- 

Computing the relevant volumes associated with the internal edges and all the vertices 
it is possible to conclude that the Regge action for this simplicial minisuperspace is 



1% S, ' 



TT — 2 arccos 



67V3 2C - 1 

ZTT -— arccos 



+ 



N, W^-3/8 

I20V2/ 

A^3 /m^/ 



6^-2 

52 



240V^V H"^ 



3/8 



(3.15) 



where Nq, Ni and N3 all refer to A4'^. 

We see that the dependence of the action on the topology of the underlying Ai^ is 
contained in the parameters Nq and A^3, (A^i — Nq+N^). The metric dependence is obviously 
contained in ^ and 5", and its dependence on the matter field in 0j and 0^. 

Note that the previous expressions are valid for a wide variety of simplicial geometries 
with very different spatial topologies (of Ms) . Different triangulations of the 3 — D spatial 
universe have different values of A^o and N3 and consequently different actions. But as long 



as they are closed connected combinatorial manifolds the previous expressions hold and so 
the functional dependence of the action on Sj and Sf, remains the same. 

So given a certain underlying simplicial complex a * Ai^ (with fixed A^^o and N-^ ) we see 
that in order to approximate the heuristic expression ( |1.1|) by a fully computable expression 

^M^; Sb, M = nNo, N,; S; 0,] = / D^D^.e-^t^'^'^^^^'^^'^-'^'l (3.16) 

we only need to specify the integration contour C, and the measure of integration D^D(f)i. 

In our simplified models the result yielded by a contour C is not very sensitive to the 
choice of measure if we stick to the usual measures, i.e., polynomials of the squared edge 
lengths. In our case we take 



dsj S_ 

2mF 2mH^ 

In the search for the correct integration contour for our simplicial minisuperspace model 
we start by reviewing some of the main results for the general continuum case. 

Unfortunately, in the case of closed cosmologies there is as yet no known explicit prescrip- 
tion for the integration contour. So usually one takes a pragmatic view, in which we look 
for contours that lead to the desired features of the wavefunction of the universe. Following 



T^, these features are: 

• It should yield a convergent path integral 

• The resulting wavefunction should predict classical spacetime in the late universe, i.e, 
oscillating behaviour when the is well approximated by the semiclassical approxi- 
mation. 

• The resulting wavefunction should obey the diffeomorphism constraints, in particular 
the Wheeler-DeWitt equation. 

It is well known that any integration contour over real metrics would yield a wavefunction 
that does not satisfy any of these basic requirements. On the other hand, an integration 
contour over complex metrics can, if wisely chosen, lead to a wavefunction that satisfies 
those requirements. 

In the simplicial framework, complex metrics arise from complex-valued squared edge 
lengths, (p.2|). The boundary squared edge length, S, has to be real and positive for obvious 



physical reasons. But the interior squared edge length, ^, can be allowed to take complex 
values. Before we attempt to choose the correct integration contour it is essential we study 
the analytical and asymptotic properties of the action as a function of the complex variable 

3.5 Analytic Study of the Action 

By allowing ^ to be complex-valued we have to contend with a multivalued Regge action, 
( p.l5| ). The deficit angle terms have an infinite number of branches, and we must deal with 
this multivaluedness, if we are to study the behaviour of steepest descents contours, (along 
which the imaginary part of the action must remain constant). Furthermore, in order to 
obtain a continuous and single-valued action, a careful analysis of all the action's branching 
points is necessary. As expected these branching points occur at the values of ^ for which 
the simplices become degenerate, i.e, where their volumes vanish. The analytic study of 
the action should then enable us to obtain integration contours that avoid such branching 
points, and along which the action is continuous and single-valued, leading to a well defined 
wavefunction of the Universe. 



The action ( 3.15 ) is trivially analytic with respect to the variables 0j, (pi, and S. But its 
dependence on the complex-valued ^ is much more complicated. So we shall investigate the 
analytic properties of I as if it was a function of ^ only, / = /[^], the other variables acting 
as parameters. 

The function has singularities at = and = 1/3, and square root branch points 
at ^ = 1/4, 1/3 and 3/8. These branch points correspond respectively to the vanishing of 
the volume of the internal 2— simplices, 3— simplices and 4— simplices. Using 



arccosz = —i log{z + v^^— T) 

we see that ^ = 1/3 is also a logarithmic branch point, near which the action behaves 
like : 

/ ~ z2V3nJ4^] log (3^- 1) (3.18) 



The multivaluedness of I[^] associated with these branch points forces us to implement 
branch cuts in order to obtain a continuous function. In general, for terms of the type 
y/z — zq we consider a branch cut (— oo, zq]. So the branch cuts associated with the terms 



— 3/8, \Ji ~ 1/4 and \Ji ~ 1/3, altogether lead to a branch cut (— c>o, 3/8]. On the other 
hand, terms of the type arccos(z) have branch points aX z = — l,+l,oo, and usually the 
associated branch cuts are chosen as (— oo, — 1] U [l,+cx3). These terms are also infinitely 
multivalued. 

The corresponding cuts for the term arccos are (1, |] U [1, 1). On the other hand, 
associated with the term arccos — ^ — we have one cut (|, fl associated with arccosu(z), 
and another (— oo, |] associated with u{z) = — 2. 

So when we consider all these branch cuts simultaneously, we see that one way to obtain 
a continuous action / as a function of ^, is to consider a total branch cut (— oo, |]. Note that 
this also takes care of the singularity at ^ = 

Although the action then becomes a continuous function of C, in the complex plane with 
a cut (— oo, |], it is still infinitely multivalued. As usual in similar cases in order to remove 
this multivaluedness we redefine the domain where the action is defined, from the complex 
plane to the Riemann surface associated with /. The infinite multivaluedness of the action 
refiects itself in I having an infinite number of branches with different values. The Riemann 
surface is composed of an infinite number of identical sheets, C— (— oo, |], one sheet for each 
branch of I. 

We define the first sheet Ci of I[^] as the sheet where the terms in arccos(2) assume 
their principal values. So the action in the first sheet will be formally equal to the original 
expression ( |3.15| ), with positive signs taken for the square root factors. Note that with the 
first sheet defined in this way, for real ^ > 3/8 the volumes and deficit angles are all real, 
leading to a real Euclidean action for ^ E [|,+oo) on the first sheet. On the other hand, 
when ^ is real and less than 1/4 in the first sheet , the volumes become pure imaginary and 
the Euclidean action becomes pure imaginary. For all other points of this first sheet the 
action is fully complex. 

Since by ( |2.2| ) we see that the simplicial metric in each 4— simplex is real iff ^ is real, 
then the simplicial geometries built out of these 4— simplices will be real when ^ is real. 
Furthermore the corresponding eigenvalues of Qij are A = {4(^ — 3/8), 1/2, 1/2, 1/2}, ||^. 
So we see that for real ^ > 3/8 we have real Euclidean signature geometries, with real 
Euclidean action, and for real ^ < 1/4, we have real Lorentzian signature geometries with 
pure imaginary Euclidean action. 

Up to now we have been concerned with what happens in the first sheet only. When we 
continue the action in ^ around one or more branch points, we will leave this first sheet and 



emerge in some other sheet of the Riemann surface. Only a few of these other sheets are 
relevant to us. 

When the action is continued in ^ once around all finite branch points — 1/4, 1/3, 3/8), 
we reach what shall be called the second sheet . It is easy to conclude that the action in this 
second sheet is just the negative of the action in the first sheet. 

I%S,cl>„<t>,]^-f%S,ct>„cf>,] (3.19) 

Once in the second sheet, if wc encircle the branch points in such a way that we cross the 
branch cut, (— oo, |], between 1/4 and 1/3, we arrive at what we shall call the third sheet. 
By doing this the terms tt — 2 arccos ^^^^ ^ and arccos ||^, both flip signs and so the action 
in this third sheet is 
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If instead of crossing the cut between 1/4 and 1/3, we transverse it between 1/3 and 3/8 
then we will end up in a different sheet. However, the asymptotic behaviour of the action is 
similar to that on the third sheet, so we will not go into the details. 

3.6 Asymptotic Behaviour of the Action 

In any investigation of the correct contour of integration it is essential to know how the 
action behaves asymptotically with respect to the variable ^, (^ — > oo), because only then 
will we be able to evaluate the contribution coming from these regions. 
In the first sheet when ^ — > oo the action behaves like 



S, (pi, 4>b] ~ — jf''^^^ S{S - Scrit)\f^ 



(3.20) 



where 



and 
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The asymptotic behaviour of P for large ^ depends on whether or not 5* is larger than the 
critical value Scrit- This value will also be crucial for the classical solutions to be obtained 
in the next section. 

First note that there are two important special cases: 

• The pure gravity model, m = 0, = 0, in which we are reduced to the kind of models 
considered in 

• The massless scalar field model, m = 0. 

In both cases 74((/)j,0b) = 0, and the asymptotic behaviour becomes independent of the 
matter field sector. In both cases as ^ — > 00 we have: 



% S, ~ ^jf2SiS - S:^f)ii (3.23) 



where 
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Note that once A(0i, 0?,) > 0, then 5"™,=° > Scrit- 

For the second sheet obviously the asymptotic behaviour of the action is just the negative 
of that in the first sheet. For the third sheet the situation is slightly different because only 
some terms in the action change sign and we have that as ^ ^ 00: 

S, 0., 0,] ~ ^^^^S{S + Slll;)^ (3.25) 

where 
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As before we consider the two sub-cases where A{(j)i, (pb) = 



J^^^-=°[e, S, 0„ 0,] ~ £^^,S{S + SiT=')f^ (3.27) 
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4 Classical Solutions 



The classical simplicial geometries are the extrema of the Regge action we have obtained 
above. In our minisuperspace model there are two degrees of freedom ^, (pi. So the Regge 
equations of motion will be: 

dl , , 

irr' (4.1) 

and 

They are to be solved for the values of ^, <^j, subject to the fixed boundary data S*, 0^. 
The classical solutions will thus be of the form ^(5*, 0fe), and <\)i{S, (pb)- The solution ^(S, (pb) 
completely determines the simplicial geometry. For the general model equations( [4.1[ ), ( |4.2| ), 
take the form 

^^^[2n - a2 arccos §4] + J' ^T'^^ {(Pi - (Pbff 
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Note that we will be expressing the values of the field (p and its mass m, in Planck units 
(Z~^). Furthermore, for simplicity of notation we have introduced ai = A^^i, 02 = ^N'ijNi 
and ag = A/'3/(4V2). 

Note that these expressions were obtained using the expression for the action on the 
first sheet. Of course, since on the second sheet the action is just the negative of this, the 
equations of motion are the same. And obviously every classical solution ^j{S,(pb) located 



on the first sheet will have a counterpart ^jj of the same numerical value, but located on 
the second sheet, and so with an action of opposite sign, I[^j{S, (pb)] = —I[C,jj{S, (pb)]. So 
the classical solutions occur in pairs. 

For physical reasons the squared boundary edge length, S, has to be real and positive, 
and we will also assume that (pb is real. 

Before we go on to study the general model it is useful to investigate the two special 
cases mentioned in section 3.6, i.e., pure gravity, and massless scalar field. 

4.1 Pure Gravity 

In this case our model is reduced to a family of models already considered by ||^, where we 
have only one metric degree of freedom ^ and its associated equation of motion is 
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The condition of real positive S means that physically acceptable classical solutions occur 
only for real ^ > 3/8 and real ^ < 1/4. There are two different kinds of classical solutions 
according to whether S'^^^ is negative or positive. 

For positive S^^^ the general form of the solutions is similar to that of the S"^ x model 
in Figure 1. There are two different regimes: For every S between and Scrit, the classical 
solutions C,j{S) = C,jj{S) are real and belong to [3/8, +oo). So they correspond to Euchdean 
signature simplicial geometries and their Euclidean actions will be real (but symmetric). For 
every S > Scrit the classical solutions ^i{S) = ^jj{S) are real and belong to (— oo, 1/4]. So 
they correspond to Lorentzian signature simplicial geometries and their Euclidean actions 
are pure imaginary 

l\^j{S)] = -I[^jj{S)]=zI[^j{S)] (4.6) 

where I[^j{S)] is real. 

When Scr ■it is negative there is no Euclidean regime. All classical solutions are Lorentzian 
geometries. For each S > 0, there is a pair of solutions C,i{S) = C,ji{S) G (— oo, 1/4] on the 
first and second sheets. Their corresponding actions are pure imaginary, as in ([4.6|). 



4.2 Massless Scalar Field 



In this case, as in the general model there are two degrees of freedom ^ and 0^. The 
corresponding equations of motion are 
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and 



0. = 06 (4.8) 

If we substitute the second equation into the first we see that it reduces to the equation 
for ^ of the previous case. So the classical solutions for ^ in the massless case coincide with 
the pure gravity case, ^{S, 0j = (pb) = (.{S). So the presence of the massless scalar field does 
not change the simplicial geometries that are classical solutions. Furthermore, the actions 
associated with these classical solutions are independent of the scalar field and coincide with 
the actions associated with the previous model. 



4.3 Massive Scalar Field 

Given the classical equations of motion ( [4. 3] ) and ( [4. 41 ), introducing the second equation into 
the first yields a cubic equation on S for each value of ^, given fixed (pb- That equation is 
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2^^-1/4^- - 6^-2^ 
where K = {l/2){ml/Hf. 

The more comphcated dependence on ^ leads us to solve the classical equation in order 
to find S for each value of ^. By inverting the resulting solutions we obtain ^{S,(j)b), and 
consequently (pi{S, (j)h), via ( W^ - This being a cubic equation we expect 3 solutions for each 
value of ^. However, because of obvious physical constraints we will accept only solutions 
that are real and positive. 

Of course as in the pure gravity model, classical solutions occur in pairs, one on the first 
sheet and the other on the second sheet. They have the same numerical value but yield 
actions of opposite sign. 

Instead of presenting ^ = ^(S", 0^) as a 3— D plot it is more informative to consider the 
graphs of ^ = ^S) for several values of (pb- Figure 2 shows one such plot for a cone over the 
0^4 triangulation of the 3— sphere, in the presence of a scalar field of mass m = 1(/^^), and a 
boundary value = l(/~^)- 

Note that near ^ = 1/4 the solution is similar to a solution in the pure gravity model. 
This is to be expected because around ^ = 1/4 the curvature (Ricci scalar) term dominates 
all others, including the matter terms in the action. So for large values of S, that is, for 
C, \ , the classical solutions in the massive scalar field model will be well approximated 
by their analogues in the pure gravity model. 

Away from the singularity the situation is different. In the Euclidean regime, ^ > 3/8, 
there is always only one real positive solution of ( [4.10|) for each value of C,- Furthermore, as 
^ — » +00, S converges to a critical value dependent on the value of the mass m and boundary 
value of the scalar field (pb- 
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This way the effective critical value decreases as the values of m or 0b increase. 
In the Lorentzian regime, ^ < 1/4, the presence of matter radically changes the classical 
solutions. When ^ is sufficiently negative, i.e., ^ < ^o, there are 3 real positive solutions of 



( [4.10| ) for each value of ^. One of these branches is almost constant and rapidly converges to 
Scrit as ^ decreases. Another branch also starts at but then converges to zero as ^ — > — oo. 
Finally, the third branch, which also converges to zero as ^ ^ — oo, but continues all the 
way to ^ = 1/4, where S +00. It is this last branch of classical solutions that is relevant 
to us, because it is the only one that continues beyond Scrit all the way to S* ^ +00. 
So for S* G [0, Scrit), we will have: 

• Two pairs of real Lorentzian signature solutions 0f,) = ^7/(5', 0;,) G (—00, 1/4], 
and {S, (f)b) = 0/ ('S', 0fe) £ (—00, 1/4] with pure imaginary Euchdean actions. 

^ ^ 

• One pair of real Euchdean signature solutions C,j {S, (pb) = ^jj{S, (pb) G [3/8, +0x3), with 
real Euclidean action. 

For S > Scrit, we will have: 

• Only one pair of real solutions C,j{S,(j)b) = C,jj{S,(j)b) G {—00, 1/4] that correspond to 
Lorentzian signature simplicial metrics, and whose Euclidean actions, though symmet- 
ric, are both pure imaginary. 

i[Us,^b)] = -mn{s,^b)] = imiis,<Pb)] 

If we increase the value of m or (pb, the value of Scrit decreases to zero and eventually 
the branch associated with the Euclidean regime vanishes. The same happens to the two 
extra Lorentzian branches that are contained in [0,5*^?,*), and we are left with a situation 
similar to the one we found in the pure gravity model when S^^^ becomes < 0. That is, 
there remains only one Lorentzian branch for the entire range of S, from S = to S ^ +00. 
For other topologies like there is no Euclidean regime whatever the value of the field, 
because Scrit < S^^'^ < 0. See Figure 3. 

5 Semiclassical Approximation 

One of the main requirements on any model is that it yields a wavefunction that in the 
late universe predicts a classical (Lorentzian) spacetime, like the one we experience. Now, a 
wavefunction of the universe will predict a classical spacetime where it is well approximated 
by the semiclassical approximation associated with Lorentzian classical solutions. We shall 



see what are the conditions that lead to such a situation in the next section. For now we 
will assume that this is the case and concentrate on the semiclassical approximation. 

The semiclassical wavefunction for our model will be obtained from the full wavefunction 



-i{^,S,(pi,4>t,) 
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by assuming that the integral is dominated by the contributions of the stationary points 
of the Regge action, i.e., the classical solutions computed above. 

We shall consider that the first integration is over the complex valued ^ and then over the 
real valued (pi,. We then assume that for each pair of boundary data {S, (pb), the integral over ^ 
is dominated by the contributions coming from the classical solutions {Cfe('S', 0&)}- If these are 
real Lorentzian solutions with purely imaginary actions Ik — ^/[^/^.(jS', 0^); 0^] = iIk{S, 06, <f)i), 
then the semiclassical approximation for the integral over ^ is 
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where ' means derivative with respect to ^ and ji}. — sgn{i"). 

The remaining integrals over 0j are now Fourier-type integrals, and can be evaluated by 

the stationary phase method by the contribution coming from their stationary points which 
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are precisely the classical solutions 0^ {S, 0^) 
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where 4(5, 0b) = /[Cfc(5, 06), 0j {S, 06)]. 

When the dominant extrema are real Euclidean solutions {Cki^^ 06)}, with real Euclidean 
actions Ik{S,(f)b) then the semiclassical evaluation of the integral over ^ leads to Laplace- 
type integrals over 0j which once more are dominated by the contributions coming from the 
stationary points of /[Cjk(-S', 06), 0i], which are precisely the classical solutions 0^ {S, 06). So 
the semiclassical wavefunction will then be 
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(5.3) 



We shall be mainly interested in the behaviour of the wavefunction of the universe for 
large S, (relative to Planck's length, /) because that corresponds to the late universe we 
experience. 

The semiclassical approximations associated with pure gravity models were calculated 
in and |^ . In these models when S is large, the only classical solutions are a pair of 
real Lorentzian geometries, C,i{S) = ^ii{S) G {—oo, 1/4], with purely imaginary Euclidean 
actions, and so the semiclassical approximation associated with them yields an oscillating 
wavefunction for 5* > Scru, as desired. 

For the massless scalar field model the classical solutions have the same simplicial geom- 
etry, as in pure gravity case, since 0j = 0f,, implies that ^ki^,4>b) = ^fe('S'), where k = I, II. 
So the action associated with the classical solutions is independent of (pi, and the influence 
of the scalar field is visible only through the pre factor in (|5.3| ). 

In the general model, with a massive scalar field, we have seen that in our range of interest, 
that is, large S, there is also a pair of real Lorentzian classical solutions, for each value of S 
(and (pb). The semiclassical wavefunction associated with both solutions has an oscillating 
behaviour as required in order to predict classical spacetime for the late universe. In Figures 
4,5, and 6 we show the numerically calculated semiclassical wavefunctions associated with 
simplicial geometries that are cones over triangulations of some of the most relevant spatial 
3 — D universes, namely simple triangulations of the 3— sphere S^, of S"^ x S^, and T^, in 
the presence of massive scalar fields. The results clearly indicate an oscillatory behaviour 
for large S, signalling classical spacetime, in all cases. In the case of T^, Scrit < 0, which 
means that the semiclassical wavefunction is valid for all values of S, and not just large 
ones. Similar calculations for other triangulations of the spatial 3 — D universe yield similar 
wavefunctions, confirming the generality of the result. 

6 Steepest Descents Contours 

We will now investigate under what conditions the semiclassical approximation is a good 
approximation of the full wavefunction. Mathematically this happens when the contour of 
integration C can be distorted such that it passes as a steepest descent contour Csd through 
the classical solutions (on which we wish to base our semiclassical approximation) so that 
the integral is dominated by the contribution from the neighbourhoods of those solutions. 
In order to guarantee that this is indeed the detailed knowledge of the entire 



contour of integration is not necessary. It is sufficient to have such a knowledge near the 
extrema that supposedly dominate the integral, and prove that the integrand is sharply 
peaked around the classical solutions and that there are no other "critical points" yielding 
significant contributions to the integral. 

In our model there are two integration variables ^ and 0j and the full wavefunction of 
the universe is given by ( p.5|) . We work under the assumption that (pi is to be integrated 
over real values, and ^ over the complex Riemann surface of the Euclidean action /. In order 
to justify the validity of the semiclassical approximation for such wavefunction we need to 
prove that the integral over ^, can be calculated as a steepest descent (SD) integral for all 
the relevant values of 0j. 

Assuming that this is true, we can replace Jq^^ dS,e~^ by its semiclassical approximation 
based on the relevant classical extrema. This will give rise to Laplace type integrals in 0j, 
when the extrema have real Euclidean action, and to Fourier type integrals in 0j, when 
the extrema have pure imaginary Euclidean action. These integrals can then be shown to 
be dominated by the stationary points of the integrand which coincide with the classical 
solutions 0j, where 

Mi^i =0 

d(j)i ^^^=^^ 

This justifies the validity of the semiclassical approximation to the full wavefunction. 

In general, a SD contour associated with an extremum ends up either at oo or at a 
singular point of the integrand, or at another extremum with the same value of Jm(J). We 
have seen that in all three models considered, i.e., pure gravity, massless scalar field, and 
massive scalar field, when S is big enough the only classical solutions are a pair or real 
Lorentzian solutions (^j{S,(f)b),(f)i{S,(j)b)) and (^jj{S,(f)b),(f)i{S,(j)b)), where = < 1/4. 
They are located on the first and second sheets respectively, and so have pure imaginary 
actions of opposite sign. Given that their actions are different valued no single SD path can 
go directly from one to the other extremum. On the other hand given that 



and 

where * denotes complex conjugation, we 
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see that the SD path that passes through 



will be the complex conjugate of the SD path that passes through ^j. So the total SD 
contour will always be composed of two complex conjugate sections, each passing through 
one extremum, and this together with the real analyticity of the action guarantees that the 
resulting wavefunction is real. 

6.1 Pure Gravity Model 

In the pure gravity model the only singularity is located at ^ = 1/3, where Im{I) diverges and 
so no SD contour that passes through a classical solution can end up at such a singularity. So 
in this model the SD paths passing through both classical Lorentzian cxtrcma arc condemned 
to end up at infinity. So we only have to worry about convergence of the integral when 
C ^ oo. 

The appropriate SD contour for S > Scru 

Csd(S) = {e e it: : Im[I(S, 0] = } (6.1) 

(where R is the Ricmann sheet of the action), can be shown to be made up of two sections, 
one passing through the Lorentzian extremum in the first sheet and the other its complex 
conjugate, passing through the extremum located in the second sheet. 

For the SD path associated with starting off in the first sheet the SD contour on the 
upper half of this sheet is given asymptotically by 

^S(S - S,rit)Im{fi) = ~m (6.2) 

and we can guarantee the convergence of the integral along this part of the SD contour 
because the real part of the action is asymptotically 

Re[mS)\ ~ ^,S{S-S^^)^\ (6.3) 

The SD contour then passes through G (— oo,l/4] crossing into the second sheet of the 
action. Proceeding along the SD contour we will then transverse the branch cut (— oo,3/8], 
again. When S is sufficiently large, the crossing happens between 1/4 and 1/3, and we 
emerge in the third sheet. If S is smaller, the crossing happens between 1/3 and 3/8, and 
we emerge in a different sheet, but the asymptotic behaviour of the action is similar and so 
in terms of the convergence of the integral, the two situations are equivalent. In the third 



sheet the SD contour then proceeds to infinity in the first quadrant along the curve that is 
asymptotically defined by 

^^S{S + Sillryrn{^i) = m (6.4) 
Once again the convergence is guaranteed by the asymptotic behaviour of the action 

Re[m S)] ^ ^,S{S + S^J.lr')^\ (6.5) 

Figure 7 shows a numerical computation of such an SD contour for the case ol M.^ — T^. 

The asymptotic behaviours of the two sections are such that they meet at infinity on the 
first and third sheets, thus obtaining a closed SD contour defining the wavefunction of the 
universe. 

The SD integral over is thus well approximated by the semiclassical approximation 
associated with these two Lorcntzian cxtrcma and so it yields a wavefunction that in the late 
universe (large S) predicts classical Lorentzian spacetime that is a solution of the Einstein 
GR equations, as desired. 

6.2 Massless Scalar Field model 

In the massless scalar field model, the situation is somewhat different. There is a second 
integration on the variable 0,. So if in the previous case we had a family of SD contours 
labelled by the value of S we now have a family of SD contours labelled by two variables, 
namely S and (pi — 4>h. 

Csd{S, <Pi-<Pb) = [ieR: Im[I{S, C, <Pi - <pb)] = /[?(^, <Pb),^i - <Pb)]] (6.6) 

where R is the Riemann sheet of the action. 

We must be able to obtain the SD contour in ^ for each value of S and (pi — 0b, and 
prove that it verifies the conditions mentioned above for all relevant values of (pi and S. The 
numerical computation of these steepest descent contours, for values S > S^^i"''^^, yields 
two different kinds of contours. For large enough values of S, the SD contours for each value 
of (pi — (pb, are similar to the SD contours we have encountered in the pure gravity model. 
See Figure 8. 

The fact that the asymptotic behaviour of the action is the same in both the pure gravity 
model and the massless scalar field model guarantees the convergence of the SD contour. 



However, when the value of S is smaller the influence of the singularity at ^ = dominates 
and we obtain a different kind of SD contour. In Figure 9 we show the SD contour for a 
cone over the 3— sphere triangulation a^, passing through the classical Lorentzian solution 
located in the first sheet, ^j, when 5" = 50 and — (j)b — 5 

We see that instead of the SD contour going off to infinity in the first quadrant of the 
first sheet, like in the previous model, it is now shifted towards ^ = by the singularity that 
exists there. 

This is possible because as ^ ^ along the SD contour, the imaginary part of the action 
remains constant while the real part diverges to +oo, 

Rem,s,<p„<p,)] ^ I ^ r' (6-7) 

Thus the contribution of the singularity to the SD integral is vanishing. 

As we pass through the Lorentzian classical solution e (— oo, 1/4] we emerge into the 
second sheet. Prom then on the behaviour of the SD contour is very similar to the one 
in the pure gravity model. There is a second branch cut crossing somewhere between 1/4 
and 3/8, (according to the value of S), and the SD contour then proceeds to infinity along 
the third sheet. Again, since the asymptotic behaviour of the action in the massless field 
models and in the pure gravity models is the same, the localisation of the contours and the 
behaviour of the real part of the action will be common to both models, and that guarantees 
the convergence of the SD integral in the massless scalar field models. 

Obviously as in the previous model, the total SD contour also includes a similar SD 
contour (complex conjugate) for the other classical Lorentzian solution, located on the second 
sheet, ensuring that the resulting wavef unction is real. 

6.3 Massive Scalar Field Model 

Although, as in the massless case, we now have a different SD contour for each value of S, (pb 
and 4>i, an appropriate SD contour leading to a Lorentzian classical spacetime can always be 
found, for each value of 0j. We are mainly interested in what happens for large values of S, 
for which there is a pair of Lorentzian classical solutions {^ki'^, (f>b), 4'i {S, (f>b)}, k — I, II, for 
every value of S and (pb- Figures 10, 11, and 12 show the results of our numerical computation 
of the SD contour passing through the classical Lorentzian solution on the first sheet, in 
the case of a cone over a^, 5^ x S^, and T^, for S — 100 , 06 = 1 and m — 1. 



They all have similar behaviours, and so we shall concentrate on the 0:4 case, in figure 
10. Proceeding upward from the extremum, ^j{S = 100,06 = 1) = 0.2114 the SD contour 
proceeds to infinity along the first quadrant of the plane, along the parabola 



^l±^ll^As(S - S^rit)Im{4l) = (6.8) 

The convergence of the integral along this part of the contour is guaranteed by the 
asymptotic behaviour of the action along the SD contour on the first sheet 

Re[l\i, S, ct>,, 0,)] ~ '^l±^^^s{S - S,ru)^J\i'\ (6.9) 

Moving downward from the extremum at — 0.2114, we immediately cross the branch 
cut and hence emerge onto the second sheet. Once again due to the alteration of the sign of 
the action one cannot proceed immediately to infinity. Instead the SD contour crosses the 
branch cut once more at ,^ = 0.2763 between 1/4 and 1/3, emerging onto the third sheet, 
where it finally proceeds to infinity in the first quadrant of this third sheet along 

^i±^^5(5 + 0/m(V^) (6.10) 
Once more convergence is a consequence of the asymptotic behaviour 

Re[f^\i^ S, 0„ 0,)] ~ ^ll^ll^AsiS + Si^!^^\ (6.11) 

As before, the SD contour that passes through the other extremum located in the equiva- 
lent position, — 0.2114, but in the second sheet, is the complex conjugate of the previous 
contour, and the full SD contour is then taken to be the union of these two sections. Thus 
we see that the semiclassical approximation based on these two Lorentzian classical solutions 
can be justified in all three cases by the existence of these SD contours. 

Note that although we only present the results for the most significant topologies the 
same computations can be performed for any other simplicial spacetime that is a cone over 
a closed, connected combinatorial 3— manifold. We only have to change the values of A^i and 
A^3 accordingly. 



7 Conclusions 



Having emphasised the need to consider more general spacetimes, other than the ones that 
are classical solutions of Einsteins' equations, i.e., manifold-based spacetimes, we then pro- 
ceeded to show that the contributions for the wavefunction of the universe coming from some 
of these non-manifold based spacetimes are still consistent with a wavefunction predicting 
classical spacetime for the late universe, in accordance with our everyday experience. In 
order to have a finite representation of these spacetimes we appealed to the simplicial rep- 
resentation of topological spaces and the Regge calculus formulation in order to specify the 
simplicial metric on those simplicial complexes. 

By considering only complexes of cone-like type (combinatorial conifolds), a natural 
minisuperspace approximation arises, with only one internal metric degree of freedom. To 
make it somewhat more realistic we also considered the presence of a massive scalar field. 

Proving the existence of convergent steepest descent contours of integration associated 
with Lorentzian classical solutions for each case, we justified the validity of the semiclassical 
approximation for the minisuperspace wavefunctions associated with our models. The result- 
ing wavefunctions do indeed exhibit an oscillatory behaviour for the late universe, consistent 
with the prediction of Lorentzian classical spacetime as we know it. 

These results add credibility to the generalization of the concept of history in the sum 
over histories formulation of quantum gravity to include more general, non-manifold based 
histories, in particular conifold based histories. Conifolds are particularly suitable for our 
objectives because although more general than manifolds they are still regular enough for 
all the usual basic concepts in GR to be easily extended to them. On the other hand 
their situation concerning algorithmic decidability is more favourable than that of manifolds 
because the set of closed connected 4— conifolds is known to be algorithmically decidable, 
which is still an open problem for 4— manifolds. 

Finally, it should be noted that of the three conditions that any contour of integration 
should verify, the SD contours we have calculated only explicitly obey two of them, namely, 
convergence and leading to a wavefunction predicting classical space-time when the universe 
is large. We have not concerned ourselves with what regards the satisfaction of the constraints 
implementing diffeomorphism invariance, the status of which is still not fully understood in 
the simplicial framework ||13|. However we expect the minisuperspace approximation to be 
a good testing ground for further study. 
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gure 1: Classical solution for the cone over the triangulation x 6 
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Figure 2: Classical solution for the cone over the triangulation with a 
scalar field of mass m = 1 and (f)b = 1 
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Figure 3: Classical solution for the cone over the triangulation with a 
scalar field of mass m = 1 and (f)b = 1 
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Figure 4: Semiclassical wavefunction for a cone over 0:4 with scalar field of 
mass m = 1 and (/jj = 1, for = 7. 
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Figure 5: Semiclassical wavefunction for a cone over x with scalar field 
of mass m = 1 and (f)b = 1, for H = 7. 
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Figure 6: Semiclassical wavefunction for a cone over with scalar field of 
mass m = 1 and (f)h = 1, for H = 7. 
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Figure 7: Steepest descents contour for a cone over with 5* = 50. Note 
that unlilce the next hgures we present both sections of the SD contour. The 
section associated with the extremum in the first sheet (at ^ = 0.2218), is 
represented by the circle-Uke points, and the other section passing through 
the second extremum located in the same position but in the second sheet 
is represented by the cross-Uke points. In both cases the contours cross the 
branch cut again at ^ = 0.2724, between the branch points 1/4 and 1/3 
before proceeding to infinity. 
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Figure 8: Steepest descents contour for a cone over q;4 Avith a massless scalar 
field taking value 0j — 06 = 0.1 and S = 100. Starting at infinity in the first 
sheet it passes through the real extremum at ^ = 0.2114, goes into the second 
sheet and crosses again the branch cut at ^ = 0.2788, between 1/4 and 1/3 
thus emerging onto the third sheet where it proceeds to infinity. Note that 
this behaviour is maintained for any larger value of S. 
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Figure 9: Steepest descents contour for a cone over a4 with a massless scalar 
field taking value — = 5 and S = 50. Moving upward from the real 
extremum in the first sheet located at = 0.1671, instead of going to infinity 
it is attracted into the singularity at ^ = 0. Moving doAvnward it goes into 
the second sheet later emerging onto the third sheet crossing the branch cut 
at ^ = 0.3052, between the branch points at 1/4 and 1/3, and proceeding to 
infinity 
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Figure 10: Steepest descents contour for a cone over a4 with scalar field of 
mass m = 1 with (ji^ = 1 and S = 100, when (pi — (pi, = 2. Starting at infinity 
in the first sheet it passes through the real extremum at ^ = 0.2114, goes into 
the second sheet and crosses again the branch cut at ,^ = 0.2763, between 
1/4 and 1/3 thus emerging onto the third sheet where it proceeds to infinity. 
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Figure 11: Steepest descents contour for a cone over x with scalar field 
of mass 171 = 1 with ^6 = 1 and S = 100, when — = 0.1. Starting at 
infinity in the first sheet it passes through the real extremum at ^ = 0.2356, 
goes into the second sheet and crosses again the branch cut at ^ = 0.2628, 
between 1/4 and 1/3 thus emerging onto the third sheet where it proceeds 
to infinity. 
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Figure 12: Steepest descents contour for a cone over with scalar field of 
mass 171 = 1 with = 1 and S = 100, when — (j)b = 0.1. Starting at 
infinity in the first sheet it passes through the real extremum at ^ = 0.2393, 
goes into the second sheet and crosses again the branch cut at ^ = 0.2598, 
between 1/4 and 1/3 thus emerging onto the third sheet where it proceeds 
to infinity. 
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